Mass formula for 2 dimensional flavorful mesons 



o 
o 

(N 
< 

in 

> 
(N 

o 

cn 

o 



X 



OSAMU Abe(^) (*) AND NOBUAKI WATANABE(^) 

(^) Laboratory of Physics 

Asahikawa Campus, Hokkaido University of Education 
9 Hokumoncho, Asahikawa 070-8621, Japan 



Preprint No.: HUEAP-15, hep-ph/0307229 



Summary. — We analytically and numerically investigate the 't Hooft equa- 
tions, the lowest order mesonic Light-Front Tamm-Dancoff equations for SU(Nc) 
and U(Nc)gauge theories, generalized to flavor non singlet mesons. We find the 
wave function can be well approximated by new basis functions and obtain an 
analytic and an empirical formulae for the mass of the lightest bound state. Its 
value is consistent with the precedent results. 

PACS ll.lO.Ef - Lagrangian and Hamiltonian approach. 

PACS 11. 10. St - Bound and unstable states; Bethe-Salpeter equations. 

PACS ll.15.Tk - Other nonperturbative techniques. 



The light front Tamm-Dancoff (LFTD) method [1, 2, 3] has been introduced 
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as an alternative tool to lattice gauge theory to investigate relativistic bound states 
nonperturbatively. In the LF coordinate, the physical vacuum is equivalent to the 
bare vacuum, since all constituents must have non-negative longitudinal momenta 
defined by = {k^ + k^) / -s/2. Because of this simple structure of the true vac- 
uum, we can avoid the serious problems which appeared in the Tamm-Dancoff (TD) 
approximation [4] in the equal time frame. Therefore, the TD approximation is 
commonly used in the context of the LF quantization. 

The techniques have been developed [5, 6, 7, 8, 9, 10] for solving LFTD equa- 
tions in several models such as the massive Schwinger model [11], which is the 
extension of the simplest (l+l)-dimensional QED2 [12]. Bergknoff [13] first applied 
LFTD approximation to the massive Schwinger model. In the most of above ref- 
erences, as they concentrated mainly on taking account of the higher Fock state 
contributions systematically in the context of LFTD approximation, they analyzed 
LFTD equations assuming that all the masses of quarks are degenerated in order 
to avoid complexities of numerical treatment. Mo and Perry [6] and Harada and 
coworkers [7] introduced basis functions to treat the massive Schwinger model in the 
context of LFTD approximation. One of the present author [10] generalized their 
basis functions. But, all the basis functions are applicable only in the case where all 
quark masses are degenerated. 

In the real world, as six quark have their inherent masses, there are many 
mesons consist of quark and anti-quark with different masses. In this short note, 
we will attempt to generalize basis functions so as to treat the masses of mesons 
consist of different flavors with different masses. We will neglect the contributions 
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from higher Fock states, then we are led to the generahzed 't Hooft-Bergknoff-EUer 
equation [13, 14] 

mi — 1 mi, — 1 



X 1 — X 



Here, parameter a specifies the model under consideration, i.e., a = for SU(A'c) 
and q; = 1 for U(A^c), P stands for the Hadamard's finite part. M is the dimension- 
less meson mass, m/ and m// are the dimensionless quark mass of flavor f and f. 

They are related to the couphng constant g and bare masses M and fhf as follows: 

_ 2nNcM^ _ 2nNcmj^ 

^ ^ (iV^ + a-l)^2' ^/ (iv2+a-l)^2- 

One of present authors (O.A.) [10] pointed out that the wave function can be 
expanded in terms of {x{l — x))^""*"-^ in case m^ = mf. Here, is the (n + 1)- 
th smallest positive solution of Eq. (6) which will be given bellow. Main interest 
in the present paper is to extend the basis functions so that we can treat mesons 
with mf 7^ mf. One may expect it is enough to extend above basis function to 
x^"+''(l — As we will see shortly, it is not the case. 

At first, according to 't Hooft [14], we put 

(3) ^x) ^x^{l-xf. 

The most singular part of the left hand side of Eq. (1) at the end point a; = e is 
given by —{m'j — l)e^~^. One of the right hand side of Eq. (1) is given by 

(4) -/37rcot(7r/3)e^-\ 

Here, we have used 

P -T ^ ^ B{a - l,h + l)F {2,1- a -h;2- a- x) 

Jo [v — xY 



-7rcot(7ra) [ax^'^l - xf - - x)^'^} 

(5) = fab{x), 

where B is the beta function and F denotes the Gauss's hypergeometric function. 
Thus, we are led to 

(6) m}-l + f37r cot(7r/3) = 0. 
Analogously, we also have at another end point x — 1 — e', 

(7) m/^ - 1 - P'n cot(7r/3') = 0. 

Here, we have used fabi^) — fbai^ — x). If we assume /3 ~ 0(m/) for small m/, we 
obtain 

(8) (3=—mf + 0{ml), 

TT 

We multiply both sides of Eq. (1) by $(a;) and integrate them, we have 

(9) M''B{1 + 2/3, 1 + 2(5') = {mj - 1)B{2(3, 1 + 2(3') + {m), - 1)5(1 + 2/3, 2(5) 

-7(/3, /3', /3, (3') + aB{l + 2/3, 1 + 2/3'). 

Here, 

Jx=0 Jy=o {y — x)'^ 

— —ira cot{7ra)B{a + c,b + d) 

+7r(a + b) cot{TTa)B{l + a + c,b + d) 

+B{-1 + a, 1 + b)B{l + c,l + d)x 

3^2(2, 1 - a - 6, 1 + c; 2 - a, 2 + c + d; 1). 



In the above equation, 3F2 denotes the generahzed hypergeometric function. Thus, 
for small nif and small nifi, we are led to 



= a 



TT 



As a next step, we consider a higher order correction to Eq. (3). The most 
general form of the wave function is given by 



'121 



N N-n N N-n 
ni=0 ji=0 712=0 J2=0 



where Pn and (3!^ are (n + l)-th smallest positive solution of Eq. (6) and Eq. (7), 
respectively. The reason why we cannot introduce a term other than one in Eq. (12) 
will be presented later. 

If we substitute Eq. (12) into Eq. (1), we have, at the end point x — e, 



(13) = - ^ (mj - 1 + 7r/?„, cot7r/3„J C^X^^-^'" 



+ E 

m, J=0 



,7 00 



ni n2 u\ 

fc=0n2j2=0 



-K-i)E E c^^-' V \ ' 

k=0n2j2=0 

Irr? n V" r^J-fci2 (~/^n2 ~ J2 + l)fc 

~\^r ~ ^) 1^ 2^ '-'ni n2 T] 

k=0n2,j2=0 

J + 1 CX) 

-EE ^(/^"i + ^ + 1 - ^) cot 7r/?„i V. 

ni 712 

fc=0 n2j2=0 

+ E E ^(/^^ +i2) cot7r/3.,<-t^ ^~^"-~ f + 

k=0n2,j2 



TT0 (jJ+'^-^h ( '^n2 ^2)k 

k\ 
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^ ^nin2 

k=0 ni,ji,n2,j2 

{2U1-Prn-Jl-P'n2-J2)k 

{2-l3n,-3i)kk\ 
-a5koB{l + + ji, i + P'n2+ h) 



S(/3„,+ji-l,l + /3' +J2) X 



Y(a + n) 

Here, (a)„ = ——^ is the Pochhammer symbol. The first hne in Eq. (14) 

r(a) 

vanishes automatically, because of the definition of 

Now, it becomes clear that Eq. (12) is the most general form. Suppose that we 
introduce the term like cx'^{l — x)'^' with 7 7^ /3n + j: then it requires 

(14) = cimj - 1 + 7r7 cot 7r7)e^~^ 

to hold. Thus, the coefficient c should vanish. Analogously, if 7' 7^ + j then 
c = 0. 

We also have similar equation to Eq. (14) at another end point x — 1 — e'. 
If we truncate Eq. (12) to given finite N. The total number of free parameters is 
(A'"+l)^(A'"+2)^/4. On the other hand, if we require Eq.(14) and similar equation to 
hold upto of order O(e^^-i) or O(e'^Jv-i)^ we have A'"(A'" + 3) independent equations. 
Thus, we cannot solve the equations in general. We have to reduce the degree of 
freedom. We put 

that is we assume 

N N 

(16) $(x) = ^d,x*+^'(l-x)^o+i+ ^ e„,„,x^"i(l-x)^"2, 

j=0 nl,n2=0 

with doo = 1 and cqo = 0. 

For a given m/ and m//, we put = M^. We can then solve Eq.(14) for 
dj and e^^^j in terms of Mj. We thus obtain the Mj dependent truncated wave 
function, say, ^{x;Mi). We can calculate a new mass eigenvalue Mj+i using this 
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wave function as 



(17) 




< ^Mi)\HmMi) > 
< $(M,)|$(Mi) > 



We can use Eq.(ll) as Mq . For N < 3, mass eigenvalue converges in 5 iterations. 
For N — 3 and < m/, m// < 0.1, we obtain M^'s by the use of Mathematica. The 
results are summarized in Tables I and II. We can fit them by polynomials: 



(18) M^{a = Q,m) = 1.8139(m/ + m/0 + 0.892(m/ + + 0.008m/ m// 



(19) M^{a = l,m) = 1 + 1.8092(m/ + m/') + 0.497(m/ + m//)^ + 1.564m/ m// 



We cannot proceed this procedure beyond N — 3, because we cannot calculate 
3^2(2, a, 6; c, 0?; 1) in desired precision. The results with m// = m/ are consistent 
with the previous results in [10]. 

Finally, we will discuss a so-called "2% discrepancy" problem. In U(A^c), single 
flavor model, we expect the dimensionless meson mass squared can be expanded as 
follows. 



By the use of the bosonization method. Banks and his coworker [15] found 61 — 
2exp(7£;) = 3.56214- • •. Bergknoff [13], however, found 61 = 27r/V3 = 3.62759- • •. 
Two results differ each other by 2%. Our result, Eq. (20), is consistent with 



+0.041 (m/ + m//)^ — 0.18m/ m/'(m/ + m//) + • • • 



+0.95(m/ + m//)^ — 4.40m/ m//(m/ + m//) + • • • . 



(20) 



= 1 + biuif + b2mj H . 



Bergknoff 's result rather than Banks et aUs. We may expect that bi — 27r/-\/3 
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is the minimum value in the context of 't Hooft-Bergknoff-EUer equation and higher 
Fock sector should be included to solve the problem. 
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Table I. - Numerical results for hound state mass in SU(Nc) model as a function of 
quark masses m f and m f . 



0.01 



0.02 



0.04 



0.06 



0.08 



0.10 



0.01 



0.03663 0.05522 0.09293 0.13136 0.17051 0.21036 



0.02 



0.05522 0.07398 0.11205 0.15084 0.19034 0.23056 



0.04 



0.09293 0.11205 0.15083 0.19033 0.23055 0.27149 



mf 0.06 
0.08 



0.13136 0.15084 0.19033 0.23056 0.27147 0.31312 



0.17051 



0.19034 0.23055 0.27147 0.31317 



0.35546 



0.10 



0.21036 0.23056 0.27149 0.31312 0.35546 



0.39866 



Table II. - Numerical results for bound state mass in U(Nc) model as a function of 
quark masses m / and m f . 



0.01 



0.02 



0.04 



0.06 



0.08 



0.10 



0.01 
0.02 
0.04 

ruf 0.06 
0.08 
0.10 



1.03660 1.05506 1.09233 1.13012 1.16847 1.20738 



1.05506 1.07387 1.11159 1.14988 1.18867 1.22799 

1.09233 1.11159 1.15040 1.18954 1.22922 1.26940 

1.13012 1.14988 1.18954 1.22965 1.27018 1.31122 

1.16847 1.18867 1.22922 1.27018 1.31154 1.35348 



1.20738 1.22799 



1.26940 



1.31122 



1.35348 1.39622 



